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ABSTRACT 



Circularization of late-type main-sequence binaries is usually attributed to turbulent 
convection, while that of early-type binaries is explained by resonant excitation of g 
modes. We show that the latter mechanism operates in solar-type stars also and is at 
least as effective as convection, despite inefficient damping of g modes in the radiative 
core. The maximum period at which this mechanism can circularize a binary composed 
of solar-type stars in 10 10 yr is as low as 3 days, if the modes are damped by radiative 
diffusion only and g-mode resonances are fixed; or as high as 6 days, if one allows for 
evolution of the resonances and for nonlinear damping near inner turning points. Even 
the larger theoretical period falls short of the observed transition period by a factor 
two. 



Subject headings: stars: binaries: close, spectroscopic — stars: oscillations 



1. Introduction 

Coeval spectroscopic binaries show an abrupt transition between circular orbits at short 
periods and eccentric orbits at longer periods. Furthermore it appears that P c irc increases with 
age among main-sequence binaries with roughly solar- mass components ( [Mathieu et al. 1992 ). 



These observations are explained by dissipative tidal interactions. In eccentric binaries, the 
tide felt by each star is time-dependent in a frame rotating with the fluid. This is also true when 
the stellar rotation is not synchronous with the orbit. Dissipation tends to damp the variable 
part of the tide at the expense of orbital or rotational energy, leading to circularization and 
synchronization. Since the tidal force depends strongly on the distance between the stars — and 
the tidal dissipation even more strongly — only short-period binaries synchronize and circularize. 

Zahn and collaborators have extensively analyzed two plausible mechanisms for tidal 
dissipation (cf. |Zahn 19*92 for a review and more complete references). The "Theory of the 



Equilibrium Tide" assumes instantaneous hydrostatic equilibrium in the tidal potential and 
attributes dissipation to turbulence in convection zones ( [Zahn 1966| , |Zahn 19771) . The theory 
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has been confirmed and roughly calibrated by comparison with binaries containing giant stars 
(Vcrbunt & Phinney 1996). An extension of the theory explains the small residual eccentricities 



of binary pulsars with evolved companions ( phinney 1992j ). The predicted transition period for 
solar-type binaries with ages ~ 10 10 yr is P c irc ~ 2.2 d flGoodman fc Oh 1997 , henceforth Paper 



I) if one allows for the fact that the largest eddies turn over more slowly than the binary orbit 
and hence do not contribute to tidal dissipation ( poldreich fc Nicholson 1977 , Goldreich fc Keely 



1977[) . Because of the latter effect, the empirical calibrations from giant stars cannot be applied 
directly to the main-sequence case, but even if the inefficiency of slow eddies is ignored, P c irc rises 
only to 6 d (Paper I). 

Notwithstanding theory, observation indicates that P c i rc ^ 11 — 12 d for local disk and halo 



dwarfs (Duquennoy & Mayor 1991, Latham et al. (1992) ). Because of the very great sensitivity 



of tidal circularization mechanisms to the orbital period, this is a large discrepancy between 
observation and theory. 

The second of Zahn's mechanisms, the "Theory of the Dynamical Tide" involves excitation of 
g modes in the radiative zones of stars [|Zahn 197q , [Zahn 1975| , |Zahn 1977fl . Since the orbital period 



is long compared to the dynamical time of the star, the g mode has a short radial wavelength, so 
that significant coupling between the mode and the tidal potential occurs only at the boundary 
between radiative and convective zones where the radial wavenumber of the g mode vanishes. 
This mechanism has been applied to circularization of early-type binaries and appears to be very 
satisfactory flGiuricin et al. 1984 , plaret fc Cunha 1997j ). These stars have convective cores and 



radiative envelopes. It is important for the theory that the thermal timescale of the envelope is 
short, so that the g modes are efficiently damped. 

Having concluded in Paper I that the Theory of the Equilibrium Tide is inadequate for 
main-sequence solar-type binaries, we decided to adapt the Theory of the Dynamical Tide to 
late-type systems, in which the cores are radiative and the envelopes convective. This appears not 
to have been done before, perhaps because the damping of g modes in radiative cores is known to 
be very slow, which appears at first sight to reduce the efficiency of the mechanism sharply (§3). 

Since our methods are similar to Zahn's, we omit many details and emphasize those aspects 
that are peculiar to the late-type case. §2 treats the excitation of the dynamical tide as if damping 
causes the g-mode resonances to overlap. In fact the resonances do not overlap unless nonlinear 
damping is efficient, but the calculation of §2 is easy and provides an upper limit to P c irc if the 
circularization mechanism depends upon the dynamical tide. §3 considers the (large) corrections 
that must be made for a finite star with discrete damped modes. Finally, §4 summarizes our results 
and discusses the prospects for resolving the apparent conflict between theory and observation 
concerning P c irc- 
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2. Tidal forcing of g modes 

In a frame corotating with the mean motion of the binary orbit, the time-dependent 
quadrupole potential exerted on star 1 by the epicyclic motion of star 2 is, to first order in the 
eccentricity (e), 

GM 2 2 /3^F 
~~a?~ r V 10 

{y 2 ,2(M)[7cos(wt - 20) - cos{ujt + 2(j))] - \/6y 2 ,ocoswt} , (2) 



<^ cxt (i,r, #,</>) = e _±rV— (1) 



where a is the semimajor axis of the binary orbit, lo ~ ^/G(Mi + M2)/a i is the orbital frequency, 
and spherical polar coordinates have been chosen with the axis perpendicular to the orbital plane. 

Synchronization is faster than circularization because the moment of inertia of the individual 
stars is small compared to that of the orbit. Therefore we assume that the stars corotate with 
the mean motion of the orbit, but we neglect coriolis forces so that we may separate variables in 
spherical polar coordinates. The relative error in the dissipation rate caused by this neglect is 
expected to be of order unity at most. Consequently the error in P c j rc is expected to be small, 
since the dissipation rate varies as a high power of the period. 

The fluid displacement in response to the potential ((2J) is conveniently divided into a 
hydrostatic "equilibrium tide" £ cq , which is calculated as though uo = 0, and a residual "dynamical 
tide" i dyn . In stably-stratified regions, 

Here <3?(r) is the potential of the unperturbed spherically-symmetric star, and 5<3? = 5& cx t + ^scif > 
where <5$ se if (r, 9, 4>) is the change in the potential of star 1 due to its distortion. Equations (J3j) are 
consequences of the following facts: 

(i) In hydrostatic equilibrium, the density and pressure are constant on equipotentials, and 
hence the entropy is also constant on equipotentials. 

(ii) In a stratified region, entropy serves as a lagrangian radial coordinate — that that is, it is 
conserved by fluid elements and differs from one mass shell to the next. 

In unstratified regions (convection zones), the actual fluid displacement does not follow 
equations @ in the limit uj — > 0. Where the entropy and molecular weight are uniform, vorticity 
is conserved. For a nonrotating star, therefore, V x £ = 0, whereas (in general) V x £ eq ^ if £ eq 
obeys equations (|3|). We regard equations @ as the definition of £ eq and require £ dyn ^ if a 
convection zone exists. 



For our purposes, we may ignore the selfgravity of both the equilibrium and the dynamical 
tides. The selfgravity of the dynamical tide is negligible in the convection zone because the mass 
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of the convection zone of a solar-type star is only a few percent of the stellar mass; and in the 
radiative core, the dynamical tide takes the form of high-order, low-frequency g-modes, whose 
natural frequencies are very accurately reproduced in the Cowling approximation because of their 
short radial wavelength. For the equilibrium tide, one can use the apsidal-motion constant (Claret 



fc Gimenez 1995) to estimate that the error in the tidal dissipation rate caused by neglecting 



^seif * s no more than 10%. 

One can show that throughout the star, £ dyn approximately satisfies the inhomogeneous linear 
equation 

&^r JD ) + l (^T ^ yn ) + W + 1) - l) &* = id + 1)^ " (4) 

In addition to the approximations discussed earlier (neglect of Coriolis forces and 53> S eif) we have 
discarded certain terms that are smaller than those displayed by factors of w 2 /u; 2 , since the tidal 
frequencies of interest are small compared with the the dynamical frequency of the star, 

_/gm 1 \ 1 / 2 

u * = • ( 5 ) 



V Rf 

Thus equation (|j) does not describe p modes, but only g modes. The righthand side of eq. @ is 
essentially the curl of the equilibrium tide and serves as a source for the dynamical tide in this 
formulation. 

The square of Brunt- Vaisala frequency, 

%T », IdP d ( pVrA 

is positive and ~ w 2 throughout most of the radiative core. In the convection zone, we take iV 2 
to vanish, although properly it has a small negative value there, and hence N 2 (r) — > at the 
boundary between the core and the envelope, which occurs at r c = O.712i?0 in the solar model 
of Bahcall fc Pinsonneault (1995)| (henceforth BP). The run of N(r) derived from their model is 



shown in Figure |T|. 

Throughout most of the core, the term (N/u) 2 on the lefthand side of eq. (Q) is very large, 
since we are concerned with tidal frequencies uj <C <-<->* ~ N . Consequently, £^ jn has a wavelength 
<C r and is accurately described by a WKBJ approximation. The radial wavenumber is 

1/2 

(7) 



+ l)(N 2 (r) ^ 



UJ 2 



Near r = r c , the wavenumber vanishes, the mode has a turning point, and ^ yn ( r ) can be 
approximated by a solution of Airy's equation, 

d 2 y 

-—r — xy w 0, where 
dx z 
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y 

r — r r 

x = — ^ (8) 



and the radial lengthscale 



A 



£(£ + 1) dN 2 



ui 2 r 2 dr 



-1/3 



0.025P d 2/3 i? . (9) 



The numerical value is derived from BP's solar model (Fig. 1), where is the orbital period 2tt/uj 
measured in days, and the subscript c denotes evaluation at r = r c . On the scale A, the righthand 
side of Eq. (|j) is of order (A 2 /r 2 ) compared to the lefthand side, so the Airy equation (||) has been 
made homogeneous for x < (i.e. r < r c ). 

In the convective envelope (r c < r < Ri), one must solve the full inhomogenous equation, 
because N 2 (r) ~ so that £^ yn varies on scales ~ r. In fact £^ yn has a constant sign except near 
the boundaries: 

£ r « at r = Ri; 

£ r & at r = r c . (10) 

The first of these conditions reflects the requirement that the surface be nearly an equipotential, 
because deviations from an equipotential have a characteristic frequency ~ l 1 / 2 ^* 3> oj. The 
second condition reflects the requirement that the solution in the convection zone match smoothly 
at r c onto the solution in the radiation zone, for which 

^ „ « r 9 -^L 

or or 

The strongest tidal components are quadrupolar (I = 2), and they have the form 



5<$> 0Jct = Real 



Qr 2 Y 2 , m (6,cP)e 



-iu)t 



where Q is a constant with dimensions of (frequency) 2 . The solution of eqs. (|4|) and (^) in the 
convection zone yields 

qr\ 

where a c is a dimensionless constant that depends mainly on the thickness of the convection zone, 
and the (9, cfr, t) dependence is implicit but understood to be the same as that of <5<J? e xt- 



Equations ( 11 ) and (||) determine the amplitude of the wave that is launched into the radiative 
core from r = r c . This wave propagates inward until it is reflected from a turning point near r = 
and returns to r = r c with a phase that depends sensitively on the tidal forcing frequency to. If 
the waves were entirely undamped in the core, their energy would increase secularly only when u 
coincided with a global g-mode eigenfrequency u;^ n ; otherwise the energy would saturate after a 
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time ~ (Acj) 1 determined by the distance Au; to the nearest resonance. Thus without damping, 
there is no secular effect on the binary eccentricity except at exact resonance. 

There is a broad regime of moderate wave damping in which energy is transferred steadily 
from the orbit to the star at a rate that is insensitive to the position of the resonances and 
independent of the actual damping rate, 7. On the one hand, if 7 <C u, then damping is too small 
to affect the tidal excitation of the ingoing wave. On the other hand, if 7 exceeds the reciprocal of 
the round-trip travel time between the turning points, 



group 



dr 



group I 



dk r 



dr, 



(12) 



then the outgoing wave returning to r c can be neglected. The frequency difference between 
neighboring resonances is Au; = 27ri~ I 4 up , so the condition 7 _1 <^ vrt group ensures that the 
resonances are broad enough to overlap. Hence this regime can be summarized as 



2irt 1 <C 7 <C oj (moderate damping). 



(13) 



Under the conditions (13), the solution (|lT)|)-(|l"l]) in the convection zone must match onto a 
solution of the homogeneous Airy equation (@) corresponding to an inward-propagating wave, 

y(x) oc Bi(x) — iAi(x), x < 0. 

The mechanical power carried by this wave is 



£ = _r 2 (l/3) W + i)\ 



•4/3 w ll/3 



5 


dN 2 


-1/3 




2" 


pr 


din r 




dr 





(14) 



With moderate damping, the quotient of Eq by the energy in the epicyclic motion of the relative 
orbit, 

g mi = - . MlM l eW +0(e 4 ), (15) 



p 2Mi + M 2 y J 

yields the rate of decay of the epicyclic energy. Evaluating Eq from equations ((l^) and (|i~4|), and 
adding the contributions of distinct quadrupolar spherical harmonics Yi m hi quadrature, we find 



"-circ-O 



1 de Ei 







e dt 2E t 



cpi 



2. 88cj' 





dN 2 


-I/a" 


pr 5 


din r 





5/3 



M 2 /Mi+M 2 \ 11/6 
M 2 R 2 V mT~ 



(16) 



We then specialize to the case of solar-type stars and double the circularization rate to allow 
for dissipation in both stars. In BP's solar model, dN 2 / 'dlnr ~ —50a; 2 at r = r c , and by numerical 
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integration of eqs. (Q) & (|To|), we find cr c ~ —0.64 [cf. eq. (11)]. The result for the circularization 
time is 

tdrc.0 ~ 8.0 x 10 3 P d 7 yr. (17) 

The subscript "0" serves as a reminder that this pertains to the moderate damping regime. 

In the simple quasilinear theory of this paper, the eccentricity decays exponentially but never 
entirely vanishes. For ease of reference, we define the transition period P c irc by 

^circ(-fcirc) = t/S, (18) 

where t is the age of the system, so that systems with periods < P c i rc and initial eccentricities 
&i ~ 0.3 have reached e 0.015. 



According to eqs. ( |r7| ) & (^), the circularization period of the oldest solar-type main-sequence 
stars would be -P c irc,o ~ 6.4 d in the moderately damped regime. As we show in the next section, 
however, solar-type binaries with comparable periods are not in the moderately damped regime 
unless nonlinear damping is efficient. 



3. Discrete Modes and Damping 



The derivation of the circularization rate (16) treated the g-mode excited at the edge of 



the core as though it were purely an ingoing wave. In fact, the tidally-forced disturbance is a 
superposition of discrete global modes — standing rather than traveling waves — whose frequencies 
are quantized: w mo d e £ Wi,n}, where n £ {1,2,.. .} indexes the number of radial nodes. The 
damping rate of these modes by radiative diffusion (which is the dominant linear damping 
mechanism) is small compared to the difference Aw = w^ n+ i — ui n between neighboring 



eigenfrequencies. This is quantified by the parameter a <C 1 estimated in equation (22) below. 
The damping is not in the "moderate" regime defined above — it is much weaker, at least if we 
consider only linear damping mechanisms. 

So it may fairly be asked why one should bother with the estimates of the previous section, 
rather than work with discrete modes from the outset. Part of the answer is that the tidal power 
Eq [eq. flliDl and the corresponding circularization rate i c irc,o [eq. fll6|) 1 are equal to the averages 
of the true quantities E{ui) and t c \ rc {uj) over the discrete modes, if the average is taken over a 
range ;> Aw but <C oj itself, and if the average is uniformly weighted in frequency. This is an 
example of a "sum rule." It is true because E is derived from linear theory and because if the star 
is subject to a short-lived tidal disturbance with a lifetime At <C (Aw) -1 = i gr oup/27r, then many 
discrete modes are excited, but the total dissipation is given correctly by the methods of §2 since 
the disturbance ends before the wave completes a round trip between the turning points. 

The second part of the answer is that the estimates of §2 are secure, whereas the consequences 
of the discreteness of the modes are complex and somewhat uncertain. A uniformly-weighted 
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average of t~ irc (u>) with respect to u (yielding t^ rc0 ) is not appropriate, since the former is 
extremely sensitive to the latter — varying by factors ~a~ 2 3> 1 as w varies by Auj — and since Co 
is dominated by the tidal interaction itself when u is close to a resonance. On the one hand, if 
the tidal interaction is the only cause of change in the distance min n \uo — coe, n \ from resonance, 
then one should average t c irc, not ^ rc , uniformly in u; this produces the "harmonic" mean (|27|), 
which is comparable to the maximum value of the instantaneous circularization time and is larger 
than t c irc,o by a factor a -1 . On the other hand, the distance from resonance is changed by effects 
independent of the tide. One such effect is gravitational radiation; this is negligible for orbital 
periods of interest, but not by a very large factor. More important is the change in the mode 
frequencies {w£ >n } as the star evolves and its mean density declines. Finally, there is good reason 
to expect that the damping rate should be enhanced by nonlinear processes, although these are 
difficult to estimate. 

In this section, we discuss the effects of the last paragraph sequentially. First we consider 
linear damping by radiative diffusion and derive the appropriate (harmonic) mean circularization 
rate if the only influence on the resonance relations is the tidal interaction itself. Then we allow 
for changes in the resonant frequencies by stellar evolution. Lastly, we demonstrate that the 
dynamical tide is significantly nonlinear near the center of the star, and explain why this is likely 
to cause additional damping. 



3.1. Damping by radiative diffusion 

The density perturbations of a g mode are accompanied by temperature fluctuations. 
Radiative diffusion tends to smooth out the temperature variations and damp the mode at the 
local rate 

Tidln n/dr 



7(r) « ^ X (r) 



1+ dhx{P/^)/dr\ ' ^ 

This formula assumes that uj <C N(r) so that the radial wavenumber k r S> [cf. eq. (0)]. The 
thermal diffusivity is given by 

4acT 3 ( dTdm r 5 k B y 1 T , on . 

X=o — TFT = -j ~5 <T L r, (20) 



3np 2 Cp V dr dr 2 fimn 

in which all symbols have their conventional meaning. 

The stable stratification of the solar core is due mainly to entropy gradients (dS/dr > 0) 
but is reinforced by a gradient in molecular weight (dfi/dr < 0), especially near the center. The 
bracketed factor in eq. (Il9|) reflects the fact that temperature perturbations are associated only 
with the entropy gradient. This factor is less than unity in the core of BP's solar model, but it is 
never negative (which would indicate instability to semiconvection) . 



The global damping rate is the average of the local rate weighted by the time spent by a wave 
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packet at each radius: 

dk r , . , / f dk r , , . 

— l{r)ir /j— dr . (21) 

o ' o 

We have approximated the inner and outer turning points by and r c , respectively, which is 
appropriate for u <C to* and i = 2. 



We define the attenuation (a) of the mode by 7i gr0 up [cf. eq. (|i2|)1. During a roundtrip between 
the turning points, the wave amplitude is damped by exp(— a). The attenuation determines 
whether the g-modes are effectively discrete or continuous: a/2n is the ratio between the half 
width at half maximum and the separation of neighboring g-mode resonances at a common I. 
Thus a ^ 1 is the condition for resonances to overlap. A numerical calculation in the solar model 
of BP yields 

P 



The P 4 dependence occurs because the local damping rate 7(7-) oc k 2 . oc P 2 [eq (0)], and the 
roundtrip time t gr0 up is also oc P 2 [eq. fll^)]- Clearly a <C 1 and the resonances are well separated 
at the circularization period predicted in §2 (6.4 d), and this must be taken into account when one 
estimates the average circularization rate. 

We can derive an approximate expression for the tidal dissipation rate from the requirements 

that 

(i) E(u) should reduce to Eq when averaged over a range large compared to Acj = \oJ^ n +i ~ ^e,n\i 
but small compared to uj itself. 

(ii) E(u) can be expressed as a sum of Breit-Wigner profiles centered at each resonance. 

Of course Eq itself depends on u, but only as a power-law, whereas E(uj) varies tremendously on 
the scale Au; ^ uo. Therefore the average considered in point (i) makes sense only if cu 3> Alo, or 
equivalently, if the modes with frequecies in the vicinity of uj are of large radial order n. In the 



solar model, n ~ lOOPd for I = 2 modes (Provost & Berthomieu 1986), so this condition is well 
satisfied. 

From the two conditions cited, it follows that 

where 7 and Aw vary slowly with n but can be treated as constants over ranges An <C n. In the 
same approximation, the infinite sum can be expressed in closed form as 

Eq sinh a 



E u) 



sin x + sinh (a/2) 



(24) 



(^H 11 ^-^ 1 - (25) 
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The tidal frequency is not constant during the history of a close binary. The circularization 
process itself removes the epicyclic energy (15) from the orbit, causing u> to increase at the rate 



^ ) circ ^circ(^) 

Therefore the time spent in the interval (uj,uj + duj) of orbital frequency is proportional to 
duj/E{uS), and it follows that the average of E(uS) with respect to time is equivalent to its harmonic 
average with respect to frequency: 

(E) h = I -^tE(u) I I JfL = E Unh(a). (27) 
J E(uj) I J E(lo) 

The corresponding circularization time of solar-type binaries is 

torch ~ a -1 t cirC)0 « 1.5 x 10 8 P] yr (P < 11.6 d), (28) 

and the transition period for systems 10 10 yr old is P c irc = 2.8 d. This is very much less than the 
value 6.4 d obtained in §2 by assuming that the resonances overlap. The situation is summarized 
graphically by the solid lines in Figure 0. 



3.2. Effects of stellar evolution 

Equations (^) and ( f27[) imply that the harmonic average of the dissipation rate is just twice 
dissipation rate halfway between resonances (i.e. at x = vr/2): E m \ n « (ol/2)Eq. By contrast, the 
instantaneous dissipation rate at exact resonance (x = 0) is E max ~ (2/o)Eq. To the extent that 
to oc E(uj), the time required to circularize is clearly dominated by the mininum value of \u\ due to 
the tide, which we denote cj m i n . Any other process that causes secular changes \u\ > cO m i n — or even 
a nonsecular process that slowly modulates u by more than Aw/2 — can substantially influence 
the mean circularization rate in the regime a <C 1. 

For example, gravitational radiation causes 

u) GR 1.7 X 10" yt\Mj d ' l 9j 

if the stars have equal masses and e< 1. Comparing this with equations (]26| ) and (^), one 
sees that gravitational radiation would begin to influence the mean circularization rate when 
e < 0.024P d 1/6 . 

The ratio of the instantaneous dissipation rate to its time average depends on oj through the 
combination x oc u — uj£ n . Thus changes in the resonant frequencies can be as important for (E) 
as changes in lo flSavonije Papaloizou 1983| ). During the main sequence phase, solar- type stars 



expand slightly: for example, Sackman et al. (1993) estimate that the Sun's radius has increased 
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by about 11.5% since the zero-age main sequence. If the individual resonant frequencies u>i n scale 
with the mean density of the star, then 

[itj "~3xW™ yr' (30) 



Clearly this effect contributes more strongly to x than gravitational waves, and with the same sign. 

3 /2 

It exceeds the minimum value x m i n due to tidal dissipation alone if e < 0.06P d , or equivalently if 
P > 2.9e^ d, where eo.3 = (e/0.3) ~ 1 is representative of the initial eccentricity (puquennoy &j 



Mayor 1991 ). At substantially smaller eccentricities or longer periods, the mean dissipation rate 



exceeds the harmonic value (27) by a factor « (x ev /4i; m i n ) 1//2 . This leads to 

tdrcev ~ 4 -8 x 10 9 e P d 3/2 yr. (31) 

This value of t C i rc is valid only if it is less than the harmonic mean (|28|), else evolutionary effects 
on the resonances are negligible, and greater than the moderately-damped value (|l~7"l), because 
as u;£,n — ► QOj averages over time and frequency are equivalent, and (E(u)) — > Eq. Once the 
eccentricity has been sufficiently reduced so that the prediction (^Tj) is less than (|i~7|), then the 
latter should be used. 

Since the evolutionary changes in {uj£ tTl } begin to be important only at periods comparable 
with the .Pcirc obtained from the harmonic average (~ 2.8 d), we expect that -P c irc is essentially 
unaffected, except in systems having a small initial eccentricity. But the later stages of 
circularization will proceed on the timescale ( |I7D rather than (|28|). 



3.3. Nonlinearity 

In linear theory, the turning points occur where N(r) = u> [cf. eq. ([?])]. Near r = 0, N(r) 
is approximately linear in r, and so the inner turning point r m ; n scales as P7 . Using the solar 
model of BP, we estimate that 

^ w 1.25 x 10- 3 P d _1 . (32) 
The smallness of the coefficient reflects the steep gradient of N(r) near r = (Fig. 1). 

In the WKBJ approximation, the energy flux of the wave is conserved, so that 
pr 2 N 2 Vg roU p\^ r \ 2 ~ constant between the turning points. (Actually a standing wave has 
zero net flux, since the fluxes of the ingoing and outgoing wave cancel. Nevertheless, this argument 
gives the correct scaling of the wave amplitude.) If £ r in this formula is the root-mean-square 
radial displacement averaged over time and angle at fixed radius, then the constant is Eq/4it, 



where Eq is given by eq. ( |14|) and is related to t c irc,o by equations (|iq)-(|16|). This allows us to 



estimate £>(?")• An appropriate dimensionless local measure of nonlinearity is 

k r £ r w 2.8 x l(T 4 e ( —) V2 P7 11/6 ( ^ . (33) 
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At the inner turning point (|32D, this becomes 

k r £ r (r miQ ) « 1.7 x 10 3 eP d 1/6 . (34) 

We conclude that the wave is strongly nonlinear near its inner turning point, unless e <C 10 -3 . 

Sufficiently large-amplitude waves are bound to damp. Where k r ^ r > 1, for example, the 
wave inverts the stratification during part of its cycle, so that denser (lower-entropy) material 
overlies less dense (higher-entropy) material. A Rayleigh- Taylor instability results provided that 
the inverted profile persists for more than a growth time, which translates to u < N. The latter 
condition is not satisfied at r m j n but it will be adequately satisfied at radii a few times larger, 
where the local nonlinearity is smaller than eq. (34) by a factor ~ (r m i n /r) 5 / 2 . So to be sure of this 



instability, one may take k r ^ r {r m \ n ) ^ 10 (say). If instability does occur, then the energy required 
to overturn and mix the material must come at the expense of the wave. This particular instability 
may be self-limiting, because mixing will reduce the molecular-weight gradient, hence move r m i n 
outward and reduce k r £ r (r m i n ). Subtler parametric instabilities may arise at weaker nonlinearity 
(k r £ r < 1), as discussed by [Kumar &; Goodman 1996 , which are less likely to mix the fluid. 



If these or other nonlinear damping mechanisms substantially reduce the wave reflected 
from the inner turning point, or even if they merely perturb its phase in a chaotic manner, then 
the narrow resonances pictured earlier will broaden and overlap. The dynamical tide will then 
effectively be in the moderate-damping regime (|l3|) , and in that case the more optimistic estimates 
of §2 for t c ; rc and P c j rc should apply. 



4. Discussion 

We have adapted Zahn's Theory of the Dynamical Tide to late-type stars with radiative cores 
and convective envelopes. As in the case of the Theory of the Equilibrium Tide (Paper I), we have 
two sets of estimates for the circularization time (t c irc) an d the transition period dividing circular 
from noncircular orbits (-P c irc)- 

Allowing only for the very inefficient linear damping of the tidally excited g modes by radiative 
diffusion, we find that the (harmonic mean) circularization time is given by eq. (f28|). Setting 
idrc — ► 3.3 Gyr (one third the main sequence lifetime of sunlike stars), we obtain a theoretical 
circularization period of only 

Pcirch « 2.8 ( Wr ^) d. (35) 

Although well short of the observationally inferred values (cf. §1 and below), this prediction is 
larger than the value ~ 2.2 d predicted by the Theory of the Equilibrium Tide when the latter is 
corrected for the slowness of the convective eddies relative to the tidal period (Paper I). Adding 
the two mechanisms raises the prediction only to ~ 3.3 d because both mechanisms depend so 
strongly on P. 
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Because of the combined effects of stellar evolution and nonlinear damping, P c i rc is likely to 
be larger than eq. ( |3"5| ) but still smaller than the observed values. As shown in §3.3, the linear 
theory indicates a strongly nonlinear wave amplitude near the inner turning point as long as 

1 /6 

e ;> 10~ 3 P d . If nonlinear losses substantially reduce the amplitude or randomize the phase 
of the wave reflected from the inner turning point, then we recover the "moderately-damped" 
circularization rate (|l7| ) and the corresponding transition period 

^,0 = 6.4(^-1 V7 d. (36) 

On the other hand, once e ^ 0.04 then stellar evolution dominates the rate of change of the 
distance from resonance even at exact resonance (§3.2), and one again recovers the larger 
circularization rate (|l7| ) and transition period (|36|). 

Mathieu et al. 1992 have concluded that t c i rc oc P 10 / 3 gives a good fit to the observed 
scaling of P c irc with binary age. The exponent of P in the harmonic average ( p8| ) is close to 10/3. 
However, the normalization of t c i rc that we predict is much too small to explain the large values of 
P c i rc that Duquennoy & Mayor 1991 and Latham et al. (1992) find for old stars (11 — 19 d). 



In summary, theory and observation appear to be in conflict regarding circularization of 
late-type main-sequence stars. The conflict must be resolved either by adjusting the theory, or by 
revising the conclusions drawn from the observations: 

(i) Perhaps one has overlooked some more effective mechanism of circularizing main-sequence 
binaries at long periods (^ 10 d). Tassoul and collaborators have proposed that rapid 
circularization and synchronization of binaries occurs by a mechanism analogous to Ekman 
pumping of fluids in rigid containers (Tassoul 1987, 1995; Tassoul & Tassoul 1996, and references 
therein). The reality of this mechanism is controversial on theoretical grounds ( [Rieutord 1992 , 



[Rieutord fe Zahn 1997 , [Tassoul fc Tassoul 1997| ). Empirically, Tassoul's mechanism does not 



seem to be required to explain circularization of early- type binaries ( piuricin et al. 1984 , Claret 



fe Cunha 1997) or systems containing giants ( Verbunt fc Phinney 199"6| ), but it contains a free 



parameter that Tassoul supposes to vary by orders of magnitude among spectral types. Other 
possibilities include tidally excited parametric instabilities ( Kumar fc Goodman 1996| ), and tidal 



excitation of inertial modes, whose existence depends essentially on coriolis forces (Savonije & 
Papaloizou 1997] ); but neither of these mechanisms has yet been applied to circularization of 
late-type binaries. 

(ii) Alternatively, perhaps the empirical inference that circularization occurs on the main 
sequence is incorrect. The significance of the trend of P c j rc with age is not clear, given the small 
numbers of binaries with well-determined masses, ages, and evolutionary states. Circularization 
is very sensitive to the history of the mean density of the stars, which can decrease by orders 
of magnitude in non-main-sequence phases of evolution (e.g. Verbunt fc Phinney 1996]) . There 



exist old double-lined spectroscopic systems whose members are almost certainly unevolved with 
circular orbits and periods above 10 d (cf. Paper I), but a small number of such systems might be 
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ascribed to chance. Furthermore, it is possible that circularization occurs before the main sequence 
(Zahn &; Bouchet 1989). If the trend of P c irc with age is real, it might result from a difference in 
mass or metallicity between the younger disk binaries and older halo systems that correlates with 
a difference in pre-main-sequence evolution. For example, protostellar disks may strongly affect 
orbital eccentricity (Artymowicz et al. 1991), and it is conceivable that the persistence of massive 
gaseous disks may correlate inversely with the abundance of dust and other metals in them. 

Larger samples of close, well-measured binaries might reduce these uncertainties, for example 
by allowing one to compare systems of similar mass and metallicity but different age. Even more 
helpful would be large samples of pre-main-sequence binaries, which it is not unreasonable to 
expect in the near future (pVlathieu 1994|) . 



We thank Pawan Kumar for helpful discussions. This work was supported by NASA under 
grant NAG5-2796. 
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r/R 

Fig. 1. — Solid curve: Solid curve: Brunt-Vaisala frequency N normalized to dynamical frequency 
uj* = (GMq/R 3 ) 1 / 2 vs. radius, based on solar model of BP. Dashed line: tidal frequency on same 
scale for P = 1 d, i.e. ujjuj* = 0.116. 
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Fig. 2. — Circularization time, averaged over resonances, versus binary period. Broken solid line: 
Fixed resonances and linear radiative damping, based on eqs. ([H]), (p8|), and (22). Dot-dashed 
line: extension of eq. (17) to P < 11.6 d given efficient nonlinear damping and/or relatively rapid 
evolution of resonant frequencies (§3.3, §3.2). Dashed lines: Effect of evolving resonant frequencies 
for eccentricities as marked; applicable only when between solid and dot-dashed lines [eq. ([31])]. 



